What do global p-modes tell us about banana cells? 
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Abstract. We have calculated the effects of giant convection cells also know as sectoral rolls 
or banana cells, on p-mode splitting coefficients. We use the technique of quasi-degenerate 
perturbation theory formulated by Lavely & Ritzwoller in order to estimate the frequency 
shifts. A possible way of detecting giant cells is to look for even splitting coefficients of 'nearly 
degenerate' modes in the observational data since these modes have the largest shifts. We find 
that banana cells having an azimuthal wave number of 16 and maximum vertical velocity of 180 
m s _1 cannot be ruled out from GONG data for even splitting coefficients. 



1. Introduction 

The power spectra of solar convective velocities show distinct peaks representing granules and 
supergranules but no distinct features at wavenumbers representative of mesogranules or giant 
cells (Wang 1989; Chou et al. 1991; Straus and Bonaccini 1997; Hathaway et al. 2000). Numerical 
simulations of solar convection routinely show the existence of mesogranules and giant cells 
(Miesh et al. 2008, Kapyla et al. 2010). Of particular interest is the existence of sectoral rolls 
or 'banana cells' having maximum radial velocities of 200-300 ms _1 . The giant convective cells 
have always been elusive to observations at the solar surface. In the past there have been studies 
which have failed to detect giant cell motions (LaBonte, Howard and Gilman 1981; Snodgrass 
and Howard 1984; Chiang, Petro and Fonkal 1987 ) as well as those hinting at their existence 
(Hathaway et al. 1996; Simon and Strous 1997). With the availability of Dopplergrams from 
SOHO/MDI, it became possible to study such long-lived and large scale features more reliably. 
Beck, Duvall and Scherrer (1998) were able to detect giant cells at the solar surface with large 
aspect ratio (~ 4) and velocities ~ 10 m s _1 using the MDI data. Hathaway et al. (2000) used 
spherical harmonic spectra from full disk measurements to detect long-lived power at I < 64. 

Are banana cells observed regularly in direct numerical simulations (DNS), real or are they 
an artifact of insufficient resolution? Is it possible for techniques of global helioseismology to 
throw any light on the existence of giant cells? Traditionally, only degenerate perturbation 
theory (DPT) has been used to calculate the effect of rotation on the p-modes. But the first 
order contribution from poloidal flows (of which giant cells are a particular case) calculated using 
the degenerate perturbation theory vanishes giving rise to the need to use the quasi degenerate 
perturbation theory (QDPT) which in contrast to DPT couples modes having slightly different 
unperturbed frequencies. Quasi degenerate perturbation theory was applied to calculate shifts 
due to flows and asphericity in solar acoustic frequencies by Lavely and Ritzwoller (1992). Roth 
and Stix (1999, 2003) used QDPT to calculate the effect of giant cells on p-modes and claimed 
that giant cells could be found by modeling the asymmetries and line broadening in the solar 
power spectrum. They claimed that finite line width of the multiplets would limit the detection 



of the frequency splittings to vertical velocity amplitude of 100 m s or larger. In Chatterjee 
and Antia (2009, CA09 now onwards), we used a different approach to examine if the use of 
quasi degenerate perturbation theory introduces significant corrections in the frequency shifts 
over that obtained from degenerate perturbation theory. We found the effect of rotation on the 
odd coefficients is negligible and hence using the degenerate theory is sufficient. Additionally, 
CA09 also calculated the effect of N-S asymmetric component of rotation, single cell meridional 
circulation, giant cells and banana cells on p-modes. We also showed that for higher angular 
degree of the poloidal flows, the 'nearly degenerate' modes would have very large frequency shifts 
(splitting coefficients). It is these splitting coefficients we can hope to detect in observational 
data. 

In this paper we present results for banana cells but with an angular degree twice that used 
in CA09. In Sect. 2, we present a brief description of the method for completeness. In Sect. 3 
we present the results of our calculations and draw conclusions in Sect. 4. 



2. Perturbation of p- &; f-modes by giant cells 

Any perturbation calculation first requires definition of a base state with unperturbed 
eigenfrequencies as well as eigenvectors denoted and respectively, which we shall call non- 
rotating spherically symmetric (NRSS) state. In this work we have used a standard solar model 
with the OPAL equation of state (Rogers and Nayfonov 2002) and OPAL opacities (Iglesias 
and Rogers 1996) and use the formulation due to Canuto and Mazitelli (1991) to calculate the 
convective flux. The equations of motion for a mode k with eigenfrequency ujk for a NRSS 
model and a model perturbed by addition of differential rotation and/or large scale flow can be 
respectively represented by, 

£ Sk = -pouls k , (1) 
£ s k + £is k = -p uj' k 2 s' k , (2) 

where ui' k is the perturbed frequency and = Ylk'eK a k' s k' is the perturbed eigenvector. Taking 
scalar product with sj in equation ([2]) and using the notation Tijk' = — f s^Cis^idV and the 
definition Cis^ = — 2iw re fpo(v.V)sk, we obtain the matrix eigenvalue equation, 

{ % 3k> + S k'j(^l' ~ w ref)} ak ' = ^'k ~ w ?ef) a J' ( 3 ) 

k'eK 

with eigenvalue A = (co' k 2 — k^ e f) and eigenvector Xj = {aj}. Here w re £ is a reference frequency 
which approximates oo' k . In this work we use w re f = ^>fc, the frequency of the mode being 
perturbed. For details on calculation of the matrix elements, the reader is referred to Sect. 3.2 
of CA 09. Essentially we use the Wigner-Eckart theorem (equation 5.4.1 of Edmonds 1960) which 
states that the general matrix element of any tensor perturbation operator can be expanded in 
terms of Wigner 3j symbols whose coefficients of expansion are independent of azimuthal order 
m and m 1 . 

An important thing to remember is that for two modes with frequencies uj\ and U2, the 
frequency shift is given by, 

= ^ (4) 

2vr 4vru;2A v ; 

where %i2 is the coupling matrix between the two modes and A = (W2 2 — wj). Traditionally 
in helioseismology, the frequency shift as a function of m is described in terms of splitting 
coefficients for all modes that are detected. These coefficients are defined by (e.g., Ritzwoller 
and Lavely 1991) 

Unim = co nl + J2 a [ q l) V l q {m), (5) 

q=0 



where ui n i is the mean frequency of the multiplet, r q (m) are the orthogonal polynomials of 

degree q and a^'s are the so called splitting coefficients. A multiplet with frequency u n i is 
21 + 1-fold degenerate in absence of rotation, magnetic field, poloidal flows and asphericity. 
Following earlier works (Lavely and Ritzwoller 1992; Roth, Howe and Komm 2002) we express 
the velocity field in terms of spherical harmonics. For completeness, we give the expression here 
again. 

v(r,M) = Re[ui(r)^(e,0)f + ^(r)V A Y?(0 s ^)- 

w\{r)v x V A Y?(M)]- (6) 

The quantities and u>* determine the radial profiles of the flows and is the horizontal 

gradient operator. The Re refers to using only the real part of the spherical harmonics as in, 

f[iy-*(M) + *?(M)]/2 if t is even, 
Re[Y t s {6,cp)} = I (7) 
t[Y-'(M)-Y s *(0,0)]/2 if t is odd. 

The first two terms in equation ([6]) define the poloidal component of the flow whereas the last 
term is the toroidal component. By the poloidal component, we imply the meridional and non- 
zonal toroidal flows (average over <p direction is zero) e.g., (i) the meridional circulation which 
carries mass poleward near the surface and sinks near the poles and (ii) the giant convection 
cells, respectively. These flows are also called large scale flows to distinguish them from other 
small scale flows like the turbulent eddies which are of the size smaller than the typical scale of 
global modes used in helioseismology. 

The sectoral rolls or banana cells are characterised by s = t. CA09 used s = t = 8 for 
banana cells. Visual inspection of the snapshots of vertical velocity in Fig. 1 of Miesch et al. 
(2006) as well as Fig. 2c of Kapyla et al. (2010) reveals that banana cells observed in DNS of 
stellar convection has s ~ 16, t ~ 16. In addition Miesch et al. (2008) also put the value of 
the maximum radial velocity at 200 m s . So it makes sense to repeat our calculations for a 
flow with an angular dependence Y^(9,(p). In this calculation we have used only the estimate 
of the magnitude of maximum radial velocities and their angular degree from the DNS. The 
solar structure we use for frequency shift calculations come from a spherically symmetric non- 
rotating standard solar model (NRSS model) and the rotation splittings needed are obtained 
from temporally averaged GONG data. These two quantities are very different for DNS since 
they have much shallower density stratification and sometimes give a rotation profile which is 
constant on cylinders rather than cones as for the Sun. A more realistic calculation would be 
to use the stratification and the rotation rate from the DNS to define the unperturbed state. 
Nevertheless, a first calculation using a NRSS model and GONG rotation profile can provide 
valuable insights into the possibility of such poloidal flows. 

In presence of only the poloidal flow (u>* = 0) we can apply the equation of mass conservation 
V.(po v ) = to get a relation between u*(r) and fg(r) e.g., 



r or 



p r 2 i4(r) " 

■ (8) 



Here po(r) is the density in a spherically symmetric solar model. So now it only remains to 
choose u l s {r) appropriately and v^r) will be determined by equation ([8]). We choose the radial 
profile of u s (r) as given by equation (19) of Roth, Howe & Komm (2002): 

4(fito P -r){r-r b ) 

ul(r) = { (Rto P -r b y (9) 

otherwise. 



Here r& = 0.7Rq and R top = Rq define the boundaries of region where the flow is confined. The 
coupling matrix in Eq. Q calculated for the flow defined above is proportional to uq. Hence 
using Eq. we have the frequency shift, 5u oc Uq. 

3. Results 

As explained in CA09, this calculation is little more involved as a non-zero t allows coupling 
between different rre and m! of the p-modes. According to the selection rules imposed by the 
Wigner 3j symbols, the s = 16, t = 16 flow couples the mode (n, /, m) with (re', I', m± 16). Hence 
it becomes important to take into account the effect of rotational splitting on A before calculating 
the effect of these kind of poloidal flows. Thus, the difference of the square of frequencies, A 
is no longer independent of m. The couplings within a multiplet i.e., (n,l,m) ^ (n,l,m ± 16) 
happen to be zero because of the anti-symmetry of the matrix elements (see Eqs. (16)-(18) in 
CA09). 

In the upper panel of Fig. [TJ we show the frequency shift of the mode with (n, I) = (11, 38) 
due to coupling with (10,46). From Eq. ([T]), we have the frequency shift, 5u oc 1/A. The 5v(m) 
shows a discontinuity at m = since A for the (11,38, m) T± (10, 46, m — 16) becomes zero at 
m = as shown in the lower panel. The reader is encouraged to compare this figure with Fig. 7 
of CA09 where they plot similar curves for coupling of the modes (18, 61, m) <^ (17, 69, m±8). It 
is these discontinuities which give rise to a large value of splitting coefficients, a q . Also notice the 
asymmetry of 5v{m) about m = in Fig. [TJ This means that not only the even coefficients, a% q 
but also the odd coefficients, a2 q +i are non-zero. Usual inversion procedure for rotation neglects 
giant cells and assumes that the odd splitting coefficients arise only from rotation. If there are 
additional contributions to these coefficients, rotation inversions may not give correct results. 
In CA09, we performed a 1.5d rotation inversion on odd coefficients a\ and to detect any 
discernible feature in the inverted profile. However, since the magnitude of the features in the 
inverted profile are smaller than the inversion errors we concluded that giant cells with uo < 100 
m s^ 1 do not have much effect on the rotation inversion. The splitting coefficients a\,a2 and 
CI3 as a function of lower turning point radius and normalised by corresponding errors in the 
observational data from GONG are shown in Fig. [2J One can easily compare these coefficients 
for s = 8, t = 8 and s = 16, t = 16 flows. While the Y^(6, eft) flow with an amplitude of no = 100 
m s _1 already produces a maximum 02 which is twice the observational error, a y/g 6 flow with 
an amplitude u$ = 300 m s _1 produces a maximum 02 which is half the observational error. 
A similar thing may be said about the coefficients a\ and 03 as well. However, we have not 
performed a 1.5 D inversion for the s = 16, t = 16 case. We do not see a clear signal of 'nearly 
degenerate' modes in the observational data. Hence, like in CA09, in order put an upper limit 
on the amplitude of the flows we compare the 02 values from theory with observational errors 
to calculate the confidence level of the upper limit. For the (9, <^>), the confidence limit 
estimated by maximum value of CL2/&2 is less than 0.6 for a flow with an amplitude no = 300m 
s _1 for several GONG data sets. In terms of the maximum radial velocity, ii™ ax ~ 180m s _1 . 
Thus banana cells with such amplitudes cannot be ruled out using the GONG data sets used. 
Another interesting observation from Fig. [2] is the large values of a q at r% ~ 0.7Rq for both t = 8 
and 16. We have checked that this is not due to r& = 0.7 Rq, but is an intrinsic property of the 
solar model used. The rotation profile used to calculate the rotational splittings for t 7^ giant 
cell flows may also be responsible. Nevertheless, since we also see this behaviour for s = 8, t = 
flow (see Fig. 5 of CA09), we may conclude that the modes which have their turning points near 
the base of the convection zone are coupled strongly due to such poloidal flows. 

4. Conclusions 

We have calculated the effects of giant cells with an azimuthal degree 16 on p-mode splitting 
coefficients. We use the technique outlined in Chatterjee & Antia (2009) to perform this 




Figure 1. Upper panel: The solid line gives the frequency shift because of interaction 
between modes (11,38, m) having lower turning radius rt = O.69i?0 and (10,46,m rfc 16) having 
r t = O.74i? , due to the flow with an angular dependence Yy?(9, 4>). The frequency shift 
for the mode (10,46) is given by the dashed line. Lower panel: A mim _ie/2i^ (solid line) for 
the coupling (ll,38,m) <^ (10,46,m — 16); and A m)m+ i6/2f (dashed line) for the coupling 
(ll,38,m) <^ (10,46,m + 16). These values are in uHz. 



calculation. We find that from observations we cannot rule out the existence of giant cells 
with an angular variation Yi®(9,(fi) and an amplitude of uq = 300 m s -1 since the splitting 
coefficients are less than the observational errors. It is important to remember that a q oc Uq. 
In other words we can say that a flow varying as Yg(9, (f>) in the Sun can be ruled out with a 
confidence level estimated by the maximum value of the ratio 02/02 ~ 2 which is 30 times more 
than the confidence level for a flow varying as Y^q(9,(/)) with the same amplitude no = 100 m 
s -1 . Speaking in terms of maximum vertical velocity, u r , according to Eq. (|6|), a no = 100 m 
s _1 for a s = 8, t = 8 flow corresponds to n™ ax = 50 m s" 1 whereas for a s = 16, t = 16 flow 
implies a n™ ax = 60 m s _1 . 

We do not find any evidence from global helioseismology which may point at the appearance 
of banana cells with an azimuthal degree t = 16 in the DNS to being just an artifact of 
insufficient grid resolution. But we believe that in the DNS the velocity spectra would peak 
at the wavenumbers of the banana cells since they are so visually conspicuous in contrast to 
observations where there is no clear peak at giant cell scales. This may be because we do 
not resolve very many scales smaller than the banana cells in the DNS. But it is important to 
remember that we have used very regular analytical expressions for the banana cells. In the 
Sun, the cells will not only be irregular but also have finite lifetimes. The GONG data sets are 
averaged over 108 days and this may be somewhat longer than the lifetimes of giant cells and 
hence the signal may be averaged out. We can look at shorter data sets, but in that case the 
errors would be larger. 
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Figure 2. (a) a\fo\, 0.2^2 and 03/(12 as a function of lower turning point radius rt for the 
Yg 8 (#, (/)) kind of flow with uo = 100 m s _1 . (b) Same as above but for a flow varying as Y^q [6, 4>) 
with uo = 300 m s . The &\, 02 and <r 3 are the errors in the corresponding observational 
splitting coefficients a±, CL2 and 03 from the GONG data set centered about November 2002. 
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